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ABSTRACT 

An  experimental  investigation  was  made  of  the  modulation  of  an 
acoustic  standing-wave  field  due  to  trapped,  cavitation-induced  bubbles. 
Maximum  acoustic  pressures  generated  are  on  the  order  of  0.1  bar.  It 
is  verified  that  bubbles  driven  below  resonance  are  trapped  just  above 
pressure  antinodes  of  the  standing-wave  field;  bubbles  driven  above  reso- 
nance, just  above  nodes.  Maximum  modulation  percentages  of  the  standing- 
wave  field  by  bubble  scattering  are  found  to  be  in  the  neighborhood  of  up 
to  12%.  Modulation  noise  spectrum  levels  are  recorded  and  are  found  to 
decrease  from  2.5  Hz  to  100  Hz  with  slopes  ranging  from  -3.2  to  -5.4  dB/ 
octave.  Less  than  1%  of  the  modulation  noise  energy  lies  in  the  region 
above  100  Hz. 


TABLE  OF  CONTENTS 

I.  INTRODUCTION 4 

II.  THEORY 6 

A.  THE  STANDING  WAVE  FIELD 6 

B.  BUBBLE  DYNAMICS  9 

III .  EXPERIMENTAL  PROCEDURE 13 

A.  MEASUREMENT  OF  NOISE 13 

B.  OBSERVATION  OF  BUBBLE  MOTION 13 

C.  CALCULATION  OF  Q 18 

D.  MEASUREMENT  OF  MODULATION  PERCENTAGE  19 

E.  MEASUREMENT  OF  PRESSURE  AMPLITUDE 19 

IV.  EXPERIMENTAL  RESULTS  20 

A.  ACOUSTIC  WAVE  FIELD 20 

B.  MODULATION  SPECTRUM 23 

C.  BUBBLE  DYNAMICS  29 

D.  MEASURED  BUBBLE  SIZES  30 

V.  CONCLUSIONS  31 

REFERENCES  32 

INITIAL  DISTRIBUTION  LIST  33 

FORM  DD  1473  34 


I.   INTRODUCTION 

This  paper  concerns  the  effects  of  cavitation-induced  bubbles  in  ampli- 
tude modulating  a  medium  intensity  (from  about  0.02  bar  (rms)  to  about  0.15 
bar  (rms))  standing  wave  field.  The  bubbles,  after  generation,  are  trapped 
at  the  nodes  or  anti -nodes  of  an  acoustic  standing  wave  field,  and  at  suffi- 
ciently high  driving  levels  dance  around.  Because  all  of  the  bubbles  gene- 
rated are  grouped  fairly  closely  in  size  around  that  value  which  would  be 
resonant  at  the  exciting  field  frequency,  they  act  as  good  sound  scatterers. 
The  constant  generation,  growth  in  size,  dancing  about  and  eventual  escape 
of  the  trapped  bubbles  causes  the  scattered  acoustic  energy  to  modulate 
randomly  the  exciting  carrier  signal.  The  purpose  of  this  research  was  to 
measure  the  frequency  spectra  for  this  random  modulation  signal. 

Preliminary  to  development  of  these  spectra,  however,  it  was  necessary 
to  investigate  the  nature  of  the  standing  wave  field,  which  led  directly 
into  another  interesting  and  as  yet  little  explored  area,  namely,  the  nature 
of  individual  bubble  motion  or  dynamics.  It  was  felt  that  a  description 
(largely  qualitative  due  to  equipment  and  time  limitations)  of  individual 
and  aggregate  bubble  behavior  would  be  useful,  both  as  an  area  for  study  in 
itself  and  as  providing  a  possible  clue  to  the  nature  of  the  observed  sound 
spectra. 

Thus,  the  investigation  resolved  itself  into  three  broad  areas,  namely, 
the  nature  of  the  standing  wave  field,  the  development  of  the  scattered 
noise  modulation  spectra,  and  the  nature  and  characteristics  of  bubble  motion 
within  the  field. 


Although  quite  a  bit  of  research  had  been  conducted  on  cavitation 
bubbles  and  cavitation  noise,  the  above  areas  have  been  relatively 
unexplored.  In  some  early  observations,  Boyle  (Pef.  1)  and  Hopwood 
(Ref.  2)  demonstrated  that  bubbles  tend  to  collect  at  the  pressure  nodes 
of  a  vertical  standing  wave  pattern.  Somewhat  later,  Blake  (Ref.  3) 
showed  that  bubbles  would  also  collect  at  pressure  anti -nodal  positions, 
provided  that  the  resonance  frequency  of  the  bubble  is  greater  than  that 
of  the  applied  sound  field.  Eller  (Ref.  4)  provided  a  theoretical 
explanation  for  the  trapping  process  and  calculated  the  minimum  pressure 
amplitude  required  to  trap  a  bubble.  Strasberg  and  Benjamin  (Ref.  5) 
observed  that  if  the  applied  field  were  strong  enough,  the  trapped 
bubbles  would  dance  about  erratically. 


II.   THEORY 

This  section  contains  a  theoretical  discussion  of  the  standing  wave 
field  and  of  certain  aspects  of  bubble  dynamics.  A  rigorous  treatment 
of  the  exceedingly  complex  modal  structure  within  the  tank  and  of  the 
minute  erratic  motions  of  the  visible  bubbles  under  the  influence  of  the 
acoustic  field  is  beyond  the  scope  of  this  paper.  Therefore,  the  following 
assumptions  are  made: 

1)  The  presence  of  pressure  release  boundaries  at  all  surfaces.  This 
is  true  at  the  water  surface,  but  only  approximately  so  at  the  sides  and 
bottom  of  the  tank,  which  are  of  thin  glass. 

2)  Perfect  perpendicularity  of  the  tank  sides,  and  that  the  water 
surface  is  parallel  to  the  tank  bottom,  neither  of  which  is  exactly  true. 

3)  The  non-existence  of  any  discrete  lines  in  the  spectrum  of  the 
noise  modulation  signal.  In  fact,  none  were  detected;  however,  weak 
frequency  components  may  have  been  present. 

A.   THE  STANDING  WAVE  FIELD 

The  following  discussion  follows  that  in  Morse  (Ref .  6) .  We  consider 
a  rectangular  tank  measuring  60  cm  by  35  cm  in  length  and  width  respec- 
tively, and  filled  to  a  depth  of  17  cm  with  water.  A  rectangular  coordinate 
system  is  established  as  in  Figure  1  below. 


The  wave  equation  in  rectangular  coordinates  is 


ifE  +  ifa  +  ^R  =  -1  ^R 
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where  x,  y,  and  z  are  spatial  coordinates,  and  p  is  acoustic  pressure. 
Solutions  to  Equation  (1)  are  of  the  form 

p  =  pasin(wjj  +  (j>x)sin(a)  I"  +  4  )sin(wz|  +  4>z)exp(-27rift)   (2) 

1         o  0   "^ 

where  f  =  ~—U  2  +   ca  2  +  co  2]  and  A  ,  4  and  4  are  determined  by  the 
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boundary  conditions.  If  we  assume  that  all  boundaries  are  pressure 

release,  the  4's  are  all  zero  and  we  have 


p  =  p  sin(w  £  )sin(u  J  )sin(o)  ^)exp  (-2-rrift)  (3) 
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where  to  =  —  m  it,  m  =1,2,3....  Similarly,  we  have  w  =  —  m  tt, 
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Each  of  these  solutions  describes  a  free  vibration  mode  of  the  enclosure, 

Thus  the  wavelength  measured  along  any  axis  must  at  all  times  be  greater 

than  -p,  as  can  be  demonstrated  bv  substituting  values  for  m  ,  m  and  m  , 
f  -  3  x   y      z 


and  calculating  the  equivalent  A  ,  A  and  A  .  The  lowest  frequency  mode 

X    y        *_> 

allowed  in  the  tank  corresponds  to  the  (111)  mode,  where  the  triplet  of 
numbers  represent  (m  m  m  )  respectively.  Applying  (4) ,  we  have  the 

j^     y  Ci 
result  that  this  mode  corresponds  to  a  frequency  of  4.18  KHz.  Below 

this  frequency  there  can  exist  no  standing  wave  structure  in  the  tank. 

At  discrete  frequencies  above  this  cutoff/  resonant  modes  can  exist, 

these  frequencies  being  determined  by  the  triplet  (m  m  m  ) .  Thus  the 

first  several  frequencies  above  the  (111)  mode,  in  ascending  order,  are 

those  corresponding  to  the  triplets  (211) ,  (311) ,  (121) ,  (221) ,  (321) , 

(131),  (411),  etc. 

At  higher  frequencies,  the  number  of  modes  that  can  exist  within  a 

given  frequency  band  becomes  greater  and  greater.  If  we  treat  the  w  , 

w  ,  a)  as  components  of  a  three  dimensional  vector  in  frequency  space  it 

can  be  shown  that  the  number  of  modes  dN  contained  within  a  frequency  band 

df  is  given  by 

dN=  ^L  (5) 

C3 

where  V  =  x  y  z  is  the  volume  of  water  in  the  tank, 
o  -^o  o  ' 

so  that  for  our  particular  tank  of  water  being  ensonified  with  a  30  KHz 
sinusoidal  signal,  we  should  have  dN  -   12  for  a  df  of  100  Hz.  Thus, 
within  a  0.3%  frequency  spread  around  30  KHz  we  should  have  about  12 
possible  modes  of  vibration. 

The  driving  transducer  forces  vibrations  in  the  tank,  with  the  result 
that  all  modes  of  vibration  can  be  excited,  but  those  modes  lying  closest 
to  the  exciting  frequency  are  more  strongly  excited  than  those  further 
from  it.  The  relative  amount  of  energy  present  in  any  mode  depends  upon 
the  Q  of  the  system,  where  Q  is  defined  by 

Q  =  j^j  (6) 

2   1 


where  f  is  the  exciting  frequency  and 

f..  and  fy   are  the  lower  and  upper  frequencies,  respectively,  at  which 
the  acoustic  pressure  amplitude  is  down  by  3  dB.  A  system  such  as  our 
tank,  with  a  relatively  high  Q  (Q  =  200-500) ,  will  have  most  of  the 
vibrational  energy  concentrated  in  only  a  few  of  the  allowed  modes 
grouped  around  the  exciting  frequency.  The  lower  the  Q  of  the  system, 
on  the  other  hand,  the  more  evenly  will  the  energy  be  distributed  over 
a  larger  number  of  allowed  modes. 

B.   BUBBLE  DYNAMICS 

A  bubble  located  in  a  standing  wave  acoustic  field  will  experience 
a  net  force  from  that  field.  This  force  is  due  to  spatial  pressure  vari- 
ations within  the  field.  If  bubble  pulsation  and  acoustic  pressure  are 
out  of  phase,  that  is,  if  the  bubble  grows  smaller  as  acoustic  pressure 
increases,  then  the  bubble  will  tend  to  be  driven  toward  a  pressure 
antinode.  If  the  reverse  condition  holds,  with  bubble  size  and  acoustic 
pressure  being  in  phase,  the  bubble  will  tend  to  be  driven  toward  a 
pressure  node.  These  phenomena  may  be  utilized  to  trap  bubbles  at 
relatively  stationary  positions  in  an  acoustic  standing  wave  field.  It 
can  be  shown  (Ref .  7)  that  the  minimum  acoustic  pressure  amplitude  at  an 
antinode  necessary  to  trap  a  bubble  in  a  standing  wave  field  in  the  z 
direction  is  given  by 

2         r(l-?2)2  +  d2,  n. 

irf       (1-C  ) 

ju         r2       poj2Rq2  and  is  equal  to  1.0  at  resonance, 

3np 
^o 

R  is  the  equilibrium  radius  of  the  bubble, 

d  is  a  damping  term  dependent  on  the  contributions  due  to  heat 

conduction,  sound  radiation,  and  viscosity, 
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n  is  a  number  between  1.0  and  1.4/  depending  on  whether  the  bubble 
oscillation  process  is  essentially  isothermal  or  adiabatic/ 
and  p  .  refers  to  the  pressure  amplitude  (Zero-to-Peak)  at  a  pressure 
antinode  of  the  wave  field. 

In  the  above  result,  we  have  assumed  that  the  wavelength  in  the  z  direction 
has  its  smallest  possible  value,  j-   .  Then  (7)  gives  us  the  minimum  pres- 
sure at  an  antinode  at  which  a  bubble  will  be  trapped  in  the  acoustic  field 
against  the  force  of  buoyancy.  Figures  2  and  3  give  ininimum  trapping  pres- 
sure vs  bubble  size  for  some  typical  frequencies.  The  significance  of  these 
figures  is  to  fix  a  rough  minimum  value  of  the  pressure  amplitudes  at  the 
anti-modes  in  the  system  during  the  experiment.  At  a  given  frequency,  the 
acoustic  pressure  necessary  to  trap  a  bubble  has  minimum  values  just  above 
and  below  resonance,  and  increases  as  we  move  away  from  resonance  in  either 
direction.  Bubbles  driven  below  resonance  will  be  trapped,  relative  to  the 
z  axis,  just  above  pressure  maxima,  while  bubbles  driven  above  resonance 
will  be  trapped  just  above  pressure  nodes. 

Regarding  the  trapping  process,  it  is  interesting  to  note  that  if  a  pure 
modal  state  or  combination  of  modal  states  could  exist  within  the  tank,  there 
could  be  no  trapping  at  frequencies  higher  than  resonance,  since  modal  posi- 
tions within  the  tank  would  be  vertical  and  horizontal  intersecting  planes, 
and  bubbles,  once  driven  to  a  vertical  plane,  would  experience  no  further 
vertical  acoustic  force,  and  would  escape.  However,  due  probably  to  the 
absorbtion  of  acoustic  energy  at  the  walls  and  various  perturbation  effects, 
such  as  the  presence  of  the  transducers  (the  driving  transducer  measured  about 
3  cm  in  diameter  and  4  cm  in  length) ,  and  non-parallelness  of  the  boundaries , 
the  theoretically  predicted  vertical  nodal  planes  were  not  detected  in  the 
sound  field,  about  which  more  will  be  said  later.  Thus,  bubbles  could  be 

trapped  above  resonance. 
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III.   EXPERIMENTAL  PROCEDURE 

The  arrangement  of  equipment  was  completely  straightforward  and  is 
illustrated  in  Figure  4.  The  sine  wave  from  the  signal  generator  was 
amplified  and  used  to  drive  a  PZT-4  transducer.  A  variable  inductor  was 
connected  in  series  to  balance  the  capacitance  of  the  transducer  and 
provide  a  more  effective  impedance  match  with  that  of  the  amplifier. 

The  acoustic  signal  was  received  on  a  small  1/8"  diameter  cylindrical 
probe  hydrophone,  given  40  dB  of  amplification,  passed  through  a  60  Hz 
filter,  then  either  viewed  directly  or  passed  through  a  rectifier  and 
detector  and  fed  to  a  1/3-octave  frequency  analyzer,  late  in  the  experi- 
ment it  was  discovered  that  a  significant  amount  of  the  carrier  was  being 
leaked  through  the  detector,  so  an  additional  low-pass  filter,  a  Kronh-Heit 
active  type,  was  installed  after  the  detector. 

The  transducers  were  installed  in  a  standard  15  gallon  fishtank, 
measuring  60  x  35  x  30  cm  in  length,  width,  and  height  respectively,  and 
filled  to  a  depth  of  17  cm  with  water.  The  driver  was  clamped  permanently 
in  one  corner  of  the  tank  and  not  touched  through  the  experiment,  while 
the  receiving  transducer  was  mounted  on  a  sliding  carriage  that  permitted 
it  to  cover  most  of  the  tank,  at  least  horizontally,  with  facility. 

A.   MEASUREMENT  OF  NOISE 

The  receiving  hydrophone  was  calibrated  at  the  start  of  the  experiment 
by  direct  comparison  with  a  calibrated  laboratory  standard.  All  pressure 
amplitudes  calculated  during  the  course  of  the  experiment  were  figured 
directly  from  the  results  of  this  calibration.  At  the  close  of  the  experi- 
ment, the  receiving  hydrophone  was  recalibrated.  The  results  of  the  two 
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calibrations  were  within  20%.  In  general,  the  pressure  measurements  con- 
tain the  largest  percentage  errors  of  any  of  the  quantities  calculated 
during  the  experiment.  Their  accuracy  (pressure  amplitude  measurements) 
is  estimated  at  ±  20%. 

A  typical  unfiltered  signal,  taken  directly  from  the  hydrophone  (after 
40  dB  of  amplification)  is  illustrated  in  Figure  5.  Our  main  interest 
lay  not  in  the  carrier  but  rather  in  the  irregular  amplitude  modulation 
of  the  carrier. 

After  the  removal  of  any  60  Hz  component  present,  the  signal  was 
passed  through  a  detector,  basically  a  low-pass  filter  and  rectifier,  with 
a  -3  dB  point  of  about  1200  Hz.  What  remains  is  the  amplitude  modulation 
signal  as  pictured  in  Figure  6. 

Very  little  could  be  conjectured  about  the  frequency  spectrum  of  the 
amplitude  modulation  signal  before  it  was  analyzed,  so  it  was  fed  into  a 
GR  119  wave  analyzer,  an  instrument  with  a  practical  low  frequency  cutoff 
of  about  20  Hz.  It  was  soon  apparent  that  by  far  the  greatest  percentage 
of  the  energy  in  the  noise  signal  lay  below  that  frequency.  Therefore, 
the  GR  119  was  replaced  with  a  GR  1564  octave  band  analyzer,  which  has  a 
low-frequency  cutoff  of  2.5  Hz,  and  the  noise  measurements  were  made  on 
that.  As  opposed  to  the  GR  119,  which  had  a  constant  width  gate  (3,10  or 
50  Hz) ,  the  GR  1564  was  capable  of  operating  in  a  1/10  or  1/3  octave  mode; 
the  latter  mode  was  used  in  the  experiment.  However,  the  GR  119  was  used 
during  several  runs  (using  the  3  Hz  gate)  to  search  for  discrete  components 
in  the  noise  signal. 

The  1/3  octave  band  level  readings  obtained  were  corrected  to  spectrum 
level  using  the  conversion: 

SPECTRUM  LEVEL  (dB)  =  BAND  LEVEL  (dB)  -  10  log  Af 
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where  Af  is  the  frequency  spread  of  the  band  in  question.  Spectrum  level 
is  then  the  average  amount  of  energy  per  hertz  bandwidth  that  one  measures 
in  any  band.  The  GR  handbook  has  a  table  of  Af  values  for  the  1/3  octave 
band  used  in  the  experiment. 

B.   OBSERVATION  OF  BUBBLE  MOTION 

The  original  experiment  was  not  set  up  for  the  purpose  of  investigating 
bubble  dynamics  within  the  tank,  and  the  apparatus  and  general  setup  did 
not  lend  themselves  to  the  taking  of  precise,  quantitative  data  in  this  area. 
However,  the  bubble  dynamics  proved  so  interesting  in  itself  that  it  was 
decided  to  include  a  section  describing,  qualitatively,  bubble  effects,  and 
some  rough  calculations  of  bubble  size  as  related  to  trapping  pressures. 

The  bubble  generation  and  trapping  process  was  observed  through  the  tank, 
side  illuminated  against  a  black  backcloth.  This  system  provided  excellent 
contrast  and  good  viewing. 

A  bubble's  size  was  determined  from  its  velocity  of  ascent  through  the 
water  by  using  Stokes1  Law.  The  assumption  was  made  that  the  bubbles' 
shapes  remained  spherical  as  they  rose,  and  no  attempt  was  made  to  apply 
empirical  corrections  to  Stokes'  Law  calculations,  since  these  corrections 
were  small  compared  to  the  error  in  determining  the  bubble's  position  and 
its  time  of  rise.  The  procedure  used  was  to  let  the  bubble  field  stabilize 
for  a  few  minutes,  then  to  select  a  bubble,  and  by  holding  a  ruler  externally 
against  the  tank,  determine  its  vertical  position.  Parallax  errors  were 
reduced  to  as  small  a  value  as  possible  by  eyeball.  The  sound  field  was 
then  rapidly  shut  off;  a  stopwatch  being  started  at  the  same  time.  The 
bubble's  time  of  ascent  to  the  surface  was  then  determined,  the  velocity 
thereby  calculated,  and  Stokes'  Law,  in  the  following  form,  employed  to  find 
bubble  radius: 
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where  v,  the  kinematic  viscosity,  =  0.0093  poise,  appropriate  for  23°C 

g  =  980  cm/sec 

y  =  bubble  velocity 

R  =  bubble  radius  in  cm 
A  plot  of  R  vs  y  is  given  in  Figure  7  below. 
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C.   CALCULATION  OF  Q 

Typical  values  for  Q  were  calculated  by  tuning  the  frequency  generator 
to  get  a  relative  peak  in  the  carrier  signal,  then  fine  tuning  it  to  locate 
the  -3  dB  points.  This  procedure  is  valid  only  as  long  as  a  single  mode  is 
present  over  the  range  of  frequencies  in  question;  that  is,  from  the  lower 
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3  dB  point  to  the  higher  3  dB  point.  Typical  values  for  this  frequency 
spread  ranged  from  50  Hz  to  250  Hz.  Refer  to  Table  1.  Our  previous  cal- 
culations, however,  indicate  that  at  30  KHz,  we  should  expect  about  12 
vibrational  modes  per  100  Hz  bandspread.  In  actuality,  although  some 
modes  were  extremely  close  together,  too  close  to  permit  a  calculation  Q, 
there  were  numerous  modes,  which,  to  the  eye,  appeared  to  be  the  only 
significant  mode  present  over  a  bandspread  of  several  hundred  hertz,  and 
thus  permitted  calculating  values  for  Q.  It  is  not  understood  why  so 
many  fo  the  theoretically  possible  model  configurations  are  absent,  or  at 
least,  immeasurable . 

D.  MEASUREMENT  OF  MODUIATION  PERCENTAGE 

Modulation  percentage,  as  defined  here,  is  the  ratio  of  the  rms  fil- 
tered modulation  signal  to  the  rms  modulated  carrier.  The  rms  filtered 
modulation  voltage  was  measured  on  the  GR-1564  octave  band  analyzer,  using 
the  all-pass  band.  During  measurements  of  modulation  percentage,  the 
hydrophone  was  located  at  a  position  in  the  sound  field  where  the  carrier 
was  a  relative  maximum. 

E.  MEASUREMENT  OF  PRESSURE  AMPLITUDE 

The  receiving  hydrophone  was  calibrated  by  comparison  with  a  laboratory 
standard.  The  receiving  signal,  regardless  of  further  processing,  was 
given  40  dB  of  amplification.  The  actual  receiving  hydrophone  sensitivity 
was  0.01  volt  (rms)  =  0.016  bar  (rms).  However,  all  scales  on  the  graphs 
included  herein  represent  actual  meter  readings,  that  is,  with  the  addition 
of  40  dB  amplification.  Therefore,  we  have  that  1.0  volt  (rms)  =  0.016 
bar  (rms) . 
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IV.  EXPERIMEOTAL  RESULTS 

A.   ACOUSTIC  WAVE  FIELD 

The  sound  field,  in  spatial  configuration,  when  probed  along  a  line 
parallel  to  one  of  the  boundaries  at  frequencies  above  the  low-frequency 
cutoff,  demonstrated  the  expected  sinusoidal  pattern.  At  relatively  low 
frequencies  the  mode  shape  came  close  to  agreeing  with  the  theoretical 
description  of  a  pure  mode.  At  higher  frequencies  however,  the  mode  shape 
appeared  to  be  the  sum  of  several  modes  of  oscillation.  Typical  field 
patterns  are  shown  in  Figure  8  and  9.  These  patterns  were  made  at  low 
driving  levels  so  that  the  effect  of  modulation  noise  was  minimal.  The 
expected  sinusoidal  shape  is  evident,  although  the  differing  heights  of 
the  various  peaks  indicate  the  presence  of  more  than  one  mode  in  the  tank. 
Note  also  that  at  no  node  does  the  acoustic  pressure  amplitude  go  to  zero, 
indicating  that  the  walls  of  the  tank  are  not  true  non-absorptive  pressure 
release  boundaries. 

An  attempt  was  made  to  construct  the  (111)  mode  within  the  tank, 
without  success.  The  lowest  obtainable  modal  configuration  was  the  (311) 
mode,  at  a  measured  5.3  KHz. 

It  is  difficult  to  give  representative  values  of  Q.  As  the  frequency 
generator  was  swept  over  a  range  of  frequencies  from,  say,  lOKHz  through 
40  KHz,  with  the  receiving  hydrophone  placed  at  an  arbitrary  spot  in  the 
tank,  there  were  observed  hundreds  of  stronger  and  weaker  resonances.  Many 
of  these  resonances  were  so  close  to  each  other  that  the  separation  between 
them  was  less  than  the  frequency  spread  necessary  to  measure  the  -3  dB 
points.  Several  resonances  at  various  frequencies  are  tabulated  in 
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Table  1  below.  The  three  values  at  each  frequency  correspond  to 
different  positions  in  the  tank.  For  a  given  mode,  then,  the  values  for 
Q  do  not  differ  by  more  than  about  20%. 

TABLE  1:  Representative  Values  of  Q 

Q 

474 
372 
429 

226 
258 
198 

197 
178 
158 

240 
254 
218 

319 
282 
269 

In  the  theoretical  calculations  it  was  determined  that  there  should 
exist  about  12  possible  modes  within  a  100  Hz  band  centered  at  30  KHz.  In 
fact,  only  two  or  three  are  observable.  Changing  the  positions  of  the 
sending  or  receiving  hydrophones  -did  not  increase  the  number  of  observed 
modes  at  a  given  frequency. 

B.  MODUIATION  SPECTRUM 

Figure  6  shows  a  typical  output  signal,  taken  directly  from  the 
receiving  hydrophone,  and  given  40  dB  amplification.  Here  V.  =70  volts 
(rms)  which  corresponds  to  the  maximum  driving  levels  used  during  the 
experiment.  Cavitation  was  severe  and  the  bubble  field  was  stabilized. 
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The  low  frequency  modulation  is  apparent.  Figure  10  shows  a  representa- 
tive unfiltered  carrier  frequency  spectrum,  with  40  dB  of  amplif i cation . 
Here  we  are  also  driving  at  a  high  level,  and  the  modulation  percentage 
is  about  10%.  The  lower  hump  in  Figure  10  is  a  measure  of  the  60  Hz 
noise  interference  in  the  output  signal  from  the  amplifier,  while  the 
upper  hump  is  the  carrier  signal  and  its  noise  modulation  products  on 
either  side. 

Measurements  of  the  modulation  percentage  indicate  that  at  the  maxi- 
mum driving  levels  used,  the  modulation  percentage  was  on  the  order  of 
8-12%.  The  filtered  modulation  signal  corresponding  to  the  unfiltered 
signal  in  Figure  6  is  shown  in  Figure  7.  The  vertical  sensitivity  of 
the  CRO  has  been  increased  to  better  delineate  the  signal.  Reference  to 
the  horizontal  sweep  time  rate  indicated  that  a  large  amount  of  the 
signal's  energy  lies  at  very  low  frequencies  (less  than  10  Hz). 

Typical  noise  modulation  spectra  are  shown  in  Figures  11,  12,  and  13. 
As  can  be  seen,  the  spectrum  level  in  every  case  dminishes  in  magnitude 
from  the  lower  to  the  higher  frequencies.  In  no  case  did  the  amount  of 
energy  in  the  modulation  spectrum  above  80  Hz  amount  to  as  much  as  1%  of 
the  total  energy  in  the  rectified  signal.  One  apparent  discrepancy  evi- 
dent in  Figures  11,  12,  and  13  is  the  fact  that  the  total  energy  of  any 
signal  measured  was  2-3  dB  greater  than  the  sum  of  the  individual  bands 
added  together.  The  only  explanation  offered  is  that  this  could  be  an 
idiosyncrasy  of  the  GR  1564  or  might  be  a  consistent  error  in  data 
recording,  namely,  a  constant  lower-than-actual  reading  of  the  dB  meter. 
This  could  amount  to  a  considerable  error  in  the  summation  of  17  bands, 
which  was  the  number  used.  During  measurement  of  these  modulation  spectra, 
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the  maximum  carrier  rms  pressure  amplitudes  were  typically  in  the  order 
of  up  to  0.04  bars.  The  percentage  modulation  was  typically  from  4%  to 
12%  at  positions  of  maximum  carrier  amplitude. 

C.   BUBBLE  DYNAMICS 

As  mentioned  earlier,  bubble  trapping  occurs  at  positions  just  above 
either  vertical  nodes  or  antinodes,  depending  on  whether  the  bubble  reso- 
nance frequency  is  lower  or  higher  than  the  acoustic  exciting  field  fre- 
quency. The  lowest  frequency  at  which  trapping  occured  was  5.8  KHz.  At 
this  frequency,  the  vertical  pressure  amplitude  profile  consisted  of  nodes 
at  the  surface  and  bottom  of  the  tank  and  an  antinode  in  the  middle.  Thus, 
only  bubbles  below  resonance  were  trapped  (just  above  an  antinode) .  At 
sufficiently  high  driving  levels,  the  bubbles  began  to  grow  (by  rectified 
diffusion)  and  finally  escaped  to  the  surface.  The  time  from  trapping  to 
final  escape,  at  this  frequency,  was  on  the  order  of  15-30  seconds  at  the 
highest  driving  levels  used. 

At  the  driving  frequencies  used  for  most  of  the  study,  that  is,  above 
about  25  KHz,  the  vertical  pressure  pattern  was  no  longer  a  single  half- 
wave,  but  several  half -waves,  thus  permitting  nodal  as  well  as  anti-nodal 
trapping.  In  fact,  the  higher  the  driving  frequency,  the  shorter  the  time 
time  period  a  given  bubble  would  stay  small  enough  to  be  trapped  in  an 
anti-nodal  position.  At  frequencies  above  30  KHz,  the  time  period  between 
observed  bubble  generation  (due  to  cavitation)  and  the  jump  from  an  anti- 
nodal  to  a  nodal  trapped  position  was  not  more  than  several  tenths  of  a 
second  and  could  not  be  quantitatively  measured  with  the  available  equipment. 
Bubbles  trapped  above  their  resonance  frequency  (that  is,  at  nodal  positions) , 
were  generally  stable,  and  depending  on  the  driving  voltage  amplitude,  would 
maintain  their  positions  for  several  minutes  at  a  time. 
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The  number  of  trapped  bubbles  in  the  tank  at  a  given  time  was  depen- 
dent upon  the  driving  frequency,  driving  voltage  amplitude,  and  the  length 
of  time  the  field  had  been  activated.  At  an  estimate,  the  number  varied 
between  15  and  50  bubbles.  Generally,  there  seemed  to  be  no  preferred 
nodal  or  antinodal  positions,  and  the  trapping  process  occured  fairly  uni- 
formly throughout  the  tank.  Quantitatively,  the  all-pass  noise  modulation 
signal  was  measured  for  several  different  frequencies  and  driving  levels, 
as  a  function  of  receiving  probe  position.  Under  no  combination  of  driving 
frequency  and  amplitude  did  the  variation  within  the  tank  exceed  2  dB. 

Given  a  trapped  bubble  at  a  fixed  frequency,  the  increase  in  its  motion 
at  higher  driving  levels  was  quite  marked.  At  lower  levels,  the  trapped 
bubbles  were  almost  stationary,  merely  undergoing  small  excursions  about  an 
equilibrium  position.  As  the  driving  voltage  was  increased,  the  motion 
became  more  violent,  following  an  ellipsoidal  path.  At  the  highest  driving 
levels  obtainable,  the  bubble  motion  amounted  to  excursions  of  as  much  as 
a  centimeter  (estimated  by  eye ) .  At  these  levels ,  the  acoustic  pressure 
amplitude  amounted  to  on  the  order  of  0.15  bar  (rms) . 

D.  MEASURED  BUBBLE  SIZES 

Bubble  sizes  calculated  tlirough  application  of  Stokes'  Law  revealed  that 
the  smallest  measured,  trapped  bubble  had  a  radius  of  85  microns,  in  an 
exciting  field  of  30.744  KHz.  Resonance  size  at  this  frequency  was  104 
microns.  Thus,  at  a  given  frequency,  there  was  a  large  variation  in  bubble 
size  on  each  side  of  resonance. 
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V.  CONCLUSIONS 

1.  Bubbles  trapped  in  a  stationary  acoustic  field  cause  scattering  of 
acoustic  energy  and  can  amplitude  modulate  the  ensonifying  acoustic  field. 

2.  Bubbles  can  be  trapped  below  and  above  resonance,  near  anti -nodes 
and  nodes  of  the  standing  wave  field,  respectively. 

3.  Carrier  signal  pressure  amplitudes  of  about  0.07  bar  (rms)  produced 
severe  bubble  motion. 

4.  At  pressure  amplitudes  ranging  up  to  about  .04  bar  (rms),  the 
amplitude  modulation  percentage  was  about  8-12%  when  bubbles  were  present, 

5.  For  the  frequencies  and  pressures  observed,  the  modulation  spectra 
levels  due  to  the  scattered  noise  decreased  with  increasing  frequency, 
with  slopes  that  ranged  from  3.2  dB  to  5.4  dB  per  octave,  over  a  range 
of  frequencies  from  2.5  Hz  to  100  Hz. 


31 


REFERENCES 


1.  Boyle,  R.  W. ,  "Ultrasonic  Stationary  Waves " ,  Nature  120,  476-477, 

(1928). 

2.  Hopv/ood,  F.  L. ,  "Ultrasonics:  Some  Properties  of  Inaudible  Sound", 

Nature  128,  748-751  (1931). 

3.  Blake,  R.  G. ,  "Bjerknes  Forces  in  Stationary  Sound  Fields", 

J.  Acoust.  Soc.  Am.  21,  551  (1949). 

4.  Eller,  A.  I.,  "Force  on  a  Bubble  in  a  Standing  Acoustic  Wave", 

J.  Acoust.  Soc.  Am.   43,  170-171  (1968). 

5.  Strasberg,  M.  and  Benjamin,  T.  B. ,  "Excitation  of  Oscillations  in 

the  Shape  of  Pulsating  Gas  Bubbles,  Experimental  Work",  J.  Acoust. 
Soc.  Am.  30,  697  (1958). 

6.  Morse,  P.  M. ,  Vibration  and  Sound,  381-429,  McGraw-Hill,  1948. 

7.  Crum,  L.  A.  and  Eller,  A.  I.,  "rIhe  lotions  of  Bubbles  in  a 

Stationary  Sound  Field",  Acoust.  Res.  Lab.  Tech  Mem.  No.  61, 
Harvard  Univ.,  1969. 


32 


INITIAL  DISTRIBUTION  LIST 

No.  Copies 

1.  Defense  Documentation  Center  2 
Cameron  Station 

Alexandria,  Virginia  22314 

2.  Library,  Code  0212  2 
Naval  Postgraduate  School 

Monterey,  California  93940 

3.  Assistant  Professor  A.  I.  Eller,  Code  61Br  3 
Department  of  Physics 

Naval  Postgraduate  School 
Monterey,  California  93940 

4.  LT  Walter  W.  Scherf ,  USNR  1 
c/o  Don  Harris 

3680  Jewell  Street 

Pacific  Beach,  California  92109 


33 


Security  Classiftca'ion 


DOCUMENT  CONTROL  DATA  -R&D 


iSecunty   das  si  fication  of  title,    body  of  abstract  and  indexing  annotation  must  be  entered  when   the  overall  report  is   classified) 


I      ORIGINATING    ACTIVITY    (Corporate   author) 

Naval  Postgraduate  School 
Monterey,  California  93940 


21.   REPORT    SECURITY     CLASSIFICATION 


Unclassified 


26.  GROUP 


3   REPORT  TITLE 


Amplitude  Modulation  of  a  Stationary  Acoustic  Field  by  Cavitation  Bubbles 

*     DESCRIPTIVE   NOTES  (Type  of  report  end, inclusive  dales) 

Master's  Thesis;  December  1971 

5     AuTHORiSifFirs/  name,   middle   initial,    last  name) 

Walter  William  Scherf 


6   REPORT  DATE 

December  1971 


la.     TOTAL    NO.    OF    PAGES 

35 


7b.    NO.    OF    RE  FS 

7 


Sa.    CONTRACT    OR    GRANT    NO 


b.    PROJEC  T    NO 


Sa.    ORIGINATOR'S    REPORT    NUMBER(S) 


9b.    OTHER    REPORT   NO(S)  (Any  other  numbers    that  may  be  assigned 
this  report) 


10      DISTRIBUTION    STATEMENT 


Approved  for  public  release;  distribution  unlimited. 


11.    SUPPLEMENTARY    NOTES 


12.    SPONSORING    MILITARY    ACTIVITY 

Naval  Postgraduate  School 
Monterey,  California  93940 


13.  ABSTRACT 


An  experimental  investigation  was  made  of  the  modulation  of  an  acoustic 
standing-wave  field  due  to  trapped,  cavitation-induced  bubbles.  Maximum  acoustics 
pressures  generated  are  on  the  order  of  0.1  bar.  It  is  verified  that  bubbles 
driven  below  resonance  are  trapped  just  above  pressure  antinodes  of  the  standing- 
wave  field;  bubbles  driven  above  resonance,  just  above  nodes.  Maximum  modulation 
percentages  of  the  standing-wave  field  by  bubble  scattering  are  found  to  be  in 
the  neighborhood  of  up  to  12%.  Modulation  noise  spectrum  levels  are  recorded  and 
are  found  to  decrease  from  2.5  Hz  to  100  Hz  with  slopes  ranging  from  -3.2  to  -5.4 
dB/octave.  Less  than  1%  of  the  modulation  noise  energy  lies  in  the  region  above 
100  Hz. 


DD,fn°or:.91473     (PAGE  M 

S/N    0101  -807-681  1 


34 


Security  Classification 


A-3140S 


Security  Classification 


KEY     WORDS 


Modulation 
Cavitation 
Standing  Wave 
Bubbles 


>,r:.,i473  <eA«) 

.101 -  607- 68  ?  1 


35 


Security  Classification 


A- 3 1 409 


b     Oui  cJ 


265  60 


Thesis 

$332 

c.l 


133701 

Scherf  . 

Amplitude  modula- 
tion of  a  stationary 
acoustic  field  by 
cavitation  bubbles^ 


Thesis 

S332 

c.l 


133700 


Scherf 

Amplitude  modula- 
tion of  a  stationary 
acoustic  field  by 
cavitation  bubbles. 


thesS332 

Amplitude  modulation  of  a  stationary  aco 


3  2768  002  00346  9 

DUDLEY  KNOX  LIBRARY 


